Transient behavior of surface plasmon polaritons scattered at a subwavelength groove 
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We present a numerical study and analytical model of the optical near-field diffracted in the 
vicinity of subwavelength grooves milled in silver surfaces. The Green's tensor approach permits 
computation of the phase and amplitude dependence of the diffracted wave as a function of the 
groove geometry, ft is shown that the field diffracted along the interface by the groove is equivalent 
to replacing the groove by an oscillating dipolar line source. An analytic expression is derived 
' from the Green's function formalism, that reproduces well the asymptotic surface plasmon polariton 

^vq , (SPP) wave as well as the transient surface wave in the near-zone close to the groove. The agreement 

between this model and the full simulation is very good, showing that the transient "near-zone" 
regime does not depend on the precise shape of the groove. Finally, it is shown that a composite 
diffractive evanescent wave model that includes the asymptotic SPP can describe the wavelength 
■ evolution in this transient near-zone. Such a semi-analytical model may be useful for the design and 

l^-} ' optimization of more elaborate photonic circuits whose behavior in large part will be controlled by 

surface waves. 
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I. INTRODUCTION 

o ; 

The analysis of light diffraction and transmission through a slit has a long history in physical optics 0, 0] . As 
discussed by Kowarz[3|, |4[ the analysis can be separated into two problems: the boundary value problem and the 
, propagation problem. The boundary value problem concerns the determination of the field immediately at the output 
plane, and interest is usually concentrated on the boundary values in the vicinity of the slit. The propagation problem 
involves determination of the field at points in the halfspace beyond the output plane in the near- and far-field. Recent 
interest in light transmission through subwavelength periodic structures with subwave length p itch j5, 6] has stimulated 
some experimental 0, H IE EE E3 and many theoretical studies 0, EE EI EE EE El El El, Efl] with the aim of 
better understanding the nature of the light field at the surface (the boundary value problem) and the influence this 
surface field on light transmission (the propagation problem). 

We report here numerical simulations of single groove and slit-groove structures using a Green's tensor method to 
solve the Maxwell field equations near the subwavelength structures on the metal / free-space interface. The simulations 
are compared with recent experimental results on single slit-groove structures [E IE liol ] and essentially confirm the 
observed amplitude and phase evolution of the scattered waves as a function of groove geometry and groove-slit 
distance. We then show that the results of the full numerical simulation can be recovered by replacing the groove 
• • . structure with a simple oscillating dipole source and again applying the Green's tensor method to find the near- and 
far-field along the surface. This oscillating dipole picture is consistent with recent charge and field distributions at 
metal-slit and metal-groove boundaries found numerically by a finite-difference time domain (FDTD) technique [13] 
and permits calculation of both the propagating and evanescent contributions to the scattered field. Finally we also 
present a simple analytic aperture- in-opaque-screen model, in the same vein as earlier models [E IE 0]j but with a 
boundary condition that posits the SPP mode at the metal/free-space interface. Comparison of the Green's tensor 
simulations to the analytic model helps to physically interpret the numeric results in terms of surface-wave modes. 
The oscillating dipole picture, however, provides deeper insight into the essential physics of surface wave generation 
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FIG. 1: Geometry of the single groove system. 



at the groove while overcoming the limitations of any fixed-boundary-condition model. 



II. NUMERICAL SIMULATIONS 



The numerical simulations are performed with the Green's tensor method [2X| . |22j . This method is very convenient 
for the study of finite-size, two-dimensional (2D) or three-dimensional (3D) objects embedded in a multilayered 
background. It relies on the resolution of the Lippmann-Schwinger equation of the electric field. 



where G(r, r') is the Green's tensor associated with the stratified background, fcg = w 2 eoA*o is the square of the 
wave propagation parameter, and Ae(r) is the "dielectric contrast," the relative permittivity difference between the 
scatterer of volume V and the adjacent layer. The Green's tensor itself is the solution of a vector wave equation with 
a point dipole source, 



where 1 is the unit tensor. The 3D Green's tensor represents the electric field at r produced by three orthogonal unit 
dipoles located at r' in a layer of dielectric constant e. An advantageous distinguishing property of the Green's tensor 
method is that only the objects of interest need be discretized. The boundary conditions at infinity are included in the 
Green's tensor of the multilayered background. In the present case, the calculation takes a very short time because 
of the small size (some tens of nanometers) of the groove and slit. Details and extensive references concerning this 
method can be found in a recent review 231. 




The system initially studied is shown in Fig.Q] It is an empty groove milled into a metallic substrate with the 
groove profile along x and extending invariantly along y. Outside the groove, the metal/free-space interface lies in 
the Oxy plane with z axis in the vertical direction. The rectangular section of the groove is characterized by width 
w and depth t. The incident free-space electromagnetic plane wave, with E-field polarized along x, impinges on the 
groove and interface at normal incidence. Because of the reduced dimensionality of the problem, all scattered waves, 
propagating and evanescent, are restricted to the Oxz plane. 

The electric field is calculated along the line z = + , x > 0. All simulations are performed at the free-space 
wavelength A = 852 nm. In order to compare these results with previously published studies [1, [I3,|2£j, the relative 
permittivity (dielectric constant) of silver is taken to be e = —33.22 + 1.17 i. The corresponding propagation length 
against absorption is 680 /im, quite long compared to Ao, because the imaginary term in the dielectric constant is 
small for silver at this wavelength. The results of the calculation are plotted in Fig.[H The principal plot in Fig. [5] 
shows the evolution of the effective index n e g = X(x)/\q as a function of the distance from the center of the groove, 
for both E x and E z components. The dotted line indicates the effective index of the SPP guided wave, 
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FIG. 2: Evolution of the E z (solid curve) and E x (dashed curve) effective index n e s(x) = A(a;)/Ao with the distance from the 
center of the groove. The groove dimensions are: t = 100 nm, w = 100 nm. Inset: evolution of the real part of E z with x, 
along the z = + line. 



for the silver/free-space interface. The variation of the effective index is sightly different for the two components close 
to the groove, but both curves converge rapidly after 3 /im. In both cases, the effective index is larger than n spp out 
to ~ 10 /im, but for greater distances, it converges to the expected n spp . The results of the simulation are consistent 
with the measurements of Ref. |8[, that reported a value of n sur f = 1-04 ± 0.01 over a distance of ~ 6 fim. They are 
also consistent with recent finite-difference time-domain (FDTD) simulations on silver surfaces [lol ] as well as similar 
measurements and simulations on gold surfaces (TTj| . The inset of Fig. [2] shows the z-component of the electric field 
diffracted by the groove along the interface. We can clearly identify two regimes. The first extends out to ~ 3 /im along 
x and is characterized by a relatively rapid decay of the amplitude. For further distances, the amplitude decreases 
much more slowly (due to absorption) and appears constant over the displayed range. This two-step evolution is 
characteristic of a transient regime within the first few micrometers from the groove. Since the incident wave is TM 
polarized (E- field parallel to x), E z belongs only to the scattered field, and does not interfere with the incident wave. 
For this reason, the mean value of the real part of E z in the total field along the interface must be zero. This is not 
the case for the x-component, as the incident field is polarized along the x direction. 
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FIG. 3: (a) Evolution of the amplitude of the 2-component of the electric field, 18 fim away from the groove ; (b) Evolution 
of the intrinsic phase ipi nt of the z-component of the electric field diffracted by the groove (color online). 



The simulations also show that the phase and the amplitude of the diffracted wave are sensitive to the groove 
dimensions, as reported in experiments [9]. In Fig.[3Ja) and (b) are plotted the amplitude and phase evolution of the 
surface wave, at large distance (18 jam). The phase <£>i n t has been determined by comparison with a cosine function 
representing the SPP guided surface wave. It corresponds to the asymptotic phase of the diffracted wave, called the 
"intrinsic phase" in Ref. @ . The evolution of the scattered wave phase and amplitude as a function of groove depth is 
typical of a resonant phenomenon. Here, the resonance concerns standing-wave modes created inside the groove. The 
w and t dimensions of the groove may be varied so as to produce a cavity that resonates when excited by an incident 
surface wave. At resonance a vertical standing mode dominates the field distribution inside the groove. Because 
of the boundary conditions, the electric field must be almost null at the bottom of the groove. Near a resonance, 
the phase of the diffracted wave varies rapidly and passes through an inflection point, while the amplitude reaches a 
maximum. As can be observed in Fig.[3l the amplitude is maximal when the phase is almost tt. The experimental 
results @ reported an intrinsic phase value of 7r/2 at the resonance groove depth. This difference of 7r/2 between the 
experiment and the simulation arises from the fact that x- and z- components of the surface wave E-field oscillate 
in quadrature. The experiment essentially measures the intrinsic phase difference in a far-field interference pattern 
between two oscillating dipoles oriented along x: one localized at the corners of a slit and the other localized at the 
groove (see inset of Fig. 0} . Thus the experiment is sensitive to the intrinsic phase of the ^-component of the surface 
E-field while the simulation calculates the intrinsic phase of the z-component. After taking this quadrature phase 
difference into account, we see that the simulations are consistent with the measurements. The groove depth for 
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FIG. 4: Comparison between the experimental results of Ref [sj and the numerical simulations performed with a groove of 
width w = 100 and depth of t = 120 nm. The silver slab is 400-nm thick, and the slit 100-nm wide. The electric-field intensity 
was computed at the exit of the slit (black cross). 

which the resonance occurs increases with the width. For a width of w = 100 nm, the simulations yield an optimal 
depth t « 110 nm. For shallow depths, some tens of nanometers, the amplitude and phase depend weakly on the 
width. At greater depths, the amplitude becomes quite sensitive to the groove width, but the phase does not change 
dramatically for widths w > 40 nm. Clearly the phase and the amplitude of the diffracted wave are very sensitive to 
the groove geometry. As the absolute groove depth is difficult to determine experimentally, the simulations results 
for ideal geometries may differ somewhat from the nominal parameters of fabricated structures. 

In the case of the slit-groove experiments reported in Q the Green's tensor simulations produced the best agreement 
with the experimental points by considering a depth of t = 120 nm, rather than the nominal experimental depth of 
100 nm. A comparison between the simulation and measurement is plotted in Fig. HI only the initial amplitude has 
been normalized to the experimental curve. Although the experimental intensity derives from a far-field interference 
fringe and the simulated field is evaluated at the output-side plane, it is legitimate to compare the two curves because 
the far- field signal is proportional to the calculated field intensity at the output-side slit exit. We note that the same 
slit-groove calculation performed with t = 100 nm is in excellent agreement with the simulation of Lalanne et al. (20l | 
using an entirely different simulation technique. 



III. FIELD SCATTERED BY A DIPOLE ALONG THE INTERFACE 



In this section we consider the 2D field radiated by a line dipole po (rather than the 3D field radiated by a point 
dipole) located just above the metal/free-space interface, as indicated in Fig. [5] This approach has been applied by 
Lalanne and Hugonin in [20I ] to study the amplitude evolution of the scattered magnetic field. The choice of placing 
the dipole just above the surface may seem arbitrary, but it is shown in the appendix that placing the dipole just 
under the interface leads to the same conclusions. The dipole is aligned parallel to the x-axis, consistent with the 
previous groove calculations of section [TTJ For the same reasons discussed there, we only calculate the expression of 
the z component of the electric field. The dipole oscillation wavelength is 852 nm. We will use the Green's tensor 
formalism to extract a simple expression for the field just above the interface. 
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FIG. 5: The source is a dipole line located along the y axis at z = , and oriented parallel to the x axis. The field is computed 
along (x, z — + ). 



The field radiated by the dipole at a point just above the interface at a distance x is given by the equation: 

k 2 

E(r, r ) = [G (r, r ) + G 5 (r, r )] p„, 

eo 

with r{x,z) = (x, + ) and ro = (0,0 + ). We denote the couple (r,ro) as (x, 0). The 3x3 tensors Go and Gs are the 
dyadic Green's functions associated with free space and the metal/free-space surface at the considered wavelength. 
Thus, the first term represents the field directly radiated by the dipole to the observation point through free space. The 
second term represents the field radiated to the observation point after reflection from the surface. The observation 
point is displaced along i, on a line just above and parallel to the surface, running through the dipole. Due to 
symmetry of the dipole radiation pattern, the z-component of the directly radiated term along the line of observation 
points is 0, and we have: 

k 2 

E z (x) = -^[G (x,0)p + G s (x,0)p ]-e z 
eo 

= M[G 5 (x,0) Po ]-e z 
eo 

= ^Gf(x,0) Po , (3) 
eo 

where e z is the unit vector of the z axis and G z s x is the zx component of the surface Green's function, i.e. the 
z-component of the Green's function produced by a dipole aligned parallel to x. 

Although there exist approximate expressions for x small compared to the wavelength (electrostatic approximation), 
these expressions are not appropriate here because the line of observation points extends far beyond a wavelength. 
The exact expression of the surface Green's function cannot be written in closed form in direct space, but we can find 
an expression susceptible to numerical evaluation by standard methods. The Green's tensor is analytically defined in 
the half Fourier space (q,z), where q is the spatial frequency parallel to the x axis. The general expression can be 
found in reference [22^ . The Gg x component is given by: 



where R is the Fresnel reflection coefficient for TM (transverse magnetic) polarization: 
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Then, Eq. ((3]) becomes: 
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In order to interpret this last equation, consider the value of the integral without the reflection coefficient: 
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FIG. 6: Comparison between the field obtained from Eq. ([5} and with the full Green's tensor numeric simulation. 



/+oo 
dqqe lqx . 
-oo 



This function I is the derivative of the function J: 

/+oo 
dqe lqx , 
-oo 



proportional to a Dirac delta function located at x = 0. So / is the derivative of a Dirac delta function, that in fact 
represents a point dipole located at x — 0: the z component of the electric field in the direction of the dipole is 
everywhere, except at two points infinitly close where it is not defined. Thus, Eq. (0 simply states that plane waves 
diffracted by the dipole in the x direction are reflected by the surface with a factor given by the Fresnel reflection 
coefficient. 

Because qR(q) is an odd function of q: 

E z(x) = / dqqR(q)sinqx. (6) 

2vre Jo 

This integral cannot be expressed in closed form but can be computed by conventional numerical techniques. In Fig. [6] 
are compared the z components of the electric field computed with the groove simulation (for t = w = 100 nm) and 
with Eq. ©. The two curves, after proper normalization, agree very well. It might appear surprising that the phase 
of the dipolar model does not need to be adjusted compared to the groove calculation, but Fig. [3] shows that the phase 
of the wave diffracted by the groove is precisely tt for this groove geometry. The overall conclusion is that it is the 
accumulation of oscillating charges at the corners of the groove, rather than details of the groove profile itself, that 
plays a key role in the global shape of the diffracted wave a few hundreds of nanometers away from the groove center. 
The field is essentially the field diffracted by a dipole placed near the surface, its structure is determined by the fact 



that the source has a broad-band spatial frequency spectrum, and that the surface supports a long-lived, guided SPP 
mode. 

An interesting point is the role of constituent propagative and evanescent modes in the creation of the surface 
wave. In particular the composite diffracted evanescent wave (CDEW) model 0], previously invoked to interpret 
similar phenomena, considers explicitly only the evanescent part. The two contributions are difficult to extract from 
full numerical simulations but can be easily carried out with the dipole approach. The scattered field of Eq. © is 
separated in its propagative and evanescent components: 



(7) 



with: 
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dqq R(q) sin qx 



dq q R(q) sin qx. 
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Figure [7] compares the real part of these two contributions. The propagative term represents a substantial fraction, 
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FIG. 7: Contribution to the field of the radiative and evanescent components to the total wave calculated from Eq.0 (color 
online) . 



its amplitude being around 50% of the total wave a few hundreds of nanometers from the dipole location. However, 
the amplitude of the propagative component decreases much faster than the evanescent term; and the wavelength of 
the propagative part is clearly longer than the wavelength of the evanescent component. The reason is that for all 
evanescent modes, including the SPP mode, q > kg, X < Xq, and n e ff > 1. For the propagative modes q < ko, X > Xq, 
n e s < 1. 

These trends appear clearly in Fig.[8j which represents the evolution with distance from the dipole of the effective 
indices of refraction for the different contributions. The index of the evanescent part is initially larger than the 
SPP index and decreases to this value asymptotically. The index of the propagative part, initially less than unity, 
approaches n = 1 from below. This is a consequence of the fact that the propagative wave is dominated by grazing 
plane waves, for which the reflection coefficient is almost equal to I, implying q ~ ko and therefore n ~ 1. When 
these two contributions are summed, the effective index follows the curve computed from the groove simulation and 
is in good agreement with measurement 0] . 



IV. ANALYTICAL MODEL 



The dipole-on-surface model can be further simplified in order to express the diffracted wave in closed form. In 
the following, we present a simple opaque-screen analytic model, similar to that of Kowarz 4] , but in which the SPP 
wave is introduced as a boundary condition on the surface. In the Kowarz model the evanescent part of the diffracted 
wave is computed assuming the presence of a slit in an infinitely thin opaque screen, and taking into account only the 
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FIG. 8: Contribution to the effective index of the radiative and evanescent components to the total wave, and comparison with 
the effective index calculated with the groove full simulation (color online). 



evanescent components of the field just above the metal. With these assumptions, the amplitude of the evanescent 
wave along the interface is given by: 



E{x) = E 
E 



where Si is the sine integral defined as: 



<?l>|fco 



sin(gw/2) ci 

q 



(8) 



Si 



i (k (x + - Si (fcoO - , \x\ > w/2, 



Si(cv) 



sini 



dt. 



The interpretation is straightforward: the slit diffracts the incident wave into a sum of evanescent waves of spatial 
frequency q whose amplitudes are weighted by the Fourier amplitudes of the slit. The Fourier spectrum of the slit is 
a cardinal sine (sine) function. Moreover, when x — > oo, the solution of Eq. [5] is correctly approximated by: 

E d 



E(x) « — — cos (kox H — 
7T x V 2 



This model however does not reproduce correctly the result of Ref . [8( . One of the reasons is that the finite conductivity 
of the screen is not included in the Kowarz approach. If we consider again a TM-polarized wave incident on the groove, 
the SPP mode of complex wave vector q spp is excited along the interface. In that case, the corresponding wave vector is 
created by the diffraction of the incident wave. The dispersion relation of the SPP guided mode, q spp = koy/e/(e + 1), 
can be retrieved by calculating the pole of the reflection coefficient of the metal/free-space interface R(q). At Ao = 852 
nm and value of e for silver, g spp = (1.0154 + ? 5.54 10 _4 )Aio. For an evanescent wave whose wave vector is near ^spp, the 
reflection coefficient can be approximated by A/(q — q S pp), where A is a constant. We can say that the incident wave 
impinging on the groove is at first diffracted with amplitude corresponding to the Fourier spectrum of the groove, 
and is then reflected along the metallic interface with a coefficient R(q). Hence, an approximate expression of the 
evanescent wave propagating along the interface is obtained by replacing Eq. ||HJ) with: 



E(x) = E Q 



-Eq 



ko sm(qw/2) 



1 



q + q spP 



c dq 



sin(qw/2) 



1 



q 



q - 9sp P 



dq. 



(9) 



Here two poles must be inserted because the SPP wave is excited in both ±x directions. This expression can be 
simplified using the fact that the width of the groove Fourier spectrum is of the order of 1/w, whereas the width of 
the "spectral line" of the Fourier SPP mode is of the order of Q(q spp ), a thousand times narrower than the Fourier 
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FIG. 9: Parametric curve representing the function i — Ei[i(l — a)t]/2n + Ei[— i(l + a)t]/2n in the complex plane. In this 
example, a — q Bpp /ko = 1.02. 



spectrum of the groove. Hence, in q— space the groove structure is essentially a constant over the width of the SPP 
response. Changing the width of the groove will only modify the amplitude of the plane waves of wave vector q « <? S pp, 
and thus the amplitude of the diffracted wave. Hence we have: 

i— ko I r oo -i 

E(x) » E / e lqx dq + E / e iqx dq. (10) 

1 + Qspp J +ko q ~ q spp 



— oo 



For x > w/2, this expression reads: 



E(x) = 2tt Eq e iq <** x K(x), (11) 



with: 



K(x) - i - Ei ~ ggppjg) , Ei (i{-k - g sp p)x) 



2vr 2tt 
The function Ei is called exponential integral and is defined by: 



oo e -t 



Ei(z) = dt. 

J—z t 

It appears that the amplitude of the SPP wave propagating along x, exp(iq spp x), is multiplied by an envelope of 
complex value K(x). This function is represented in the complex plane for typical parameters of SPP wave vector 
in Fig. [9]^ a). The low values of x correspond to the right of the curve. When x goes toward infinity, K(x) whirls 
toward zu m = i. The strong oscillation at the beginning of the curve is due to a beating between the Ei(— q spp ) and 
the Ei(+q spp ) term. As this function has a varying phase, it will affect the wavelength of the surface wave. Figure 
[TU] compares the evolution of the surface index for x— and z-components using the index computed from the previous 
formula. This analytical model predicts the same trends as the numerical and the dipole approach: the effective index 
of the wave generated near the groove is greater than n spp but decreases and converges toward the expected SPP 
value within a range of about 10 ^m. In fact, the effective index oscillates slightly around n spp at larger distances. 
There is a good qualitative agreement with the dipole model, because the evanescent waves play the main role in the 
creation of the surface wave. The fact that the reflection coefficient is replaced by a simple pole in q spp modifies the 
time evolution of the wave amplitude (not shown), but the wavelength evolution is unchanged. The effective index is 
overestimated in the first micrometers near the groove because the radiative part is not taken into account. 



V. SUMMARY 



We have studied in detail the structure of the wave diffracted by a groove or slit milled in a metallic surface, 
illuminated by a monochromatic plane wave. First, the Green's tensor method has been used to analyze the amplitude, 
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FIG. 10: Comparison between the E z , effective index and the effective index computed from the analytical model. The groove's 
dimensions are: t — 100 nm, w = 100 nm. 



phase and frequency behavior of the surface wave in the vicinity of the groove. In this zone, the surface wave has 
a transient regime characterized by a rapid variation of the amplitude within the first 2 to 3 micrometers, and an 
increase of the surface wavelength up to the value of the SPP-wavelength in the first 15 micrometers. The phase 
and the amplitude of the scattered wave depends strongly on the groove geometry, as the incident wave excites an 
"organ-pipe" mode inside the groove. Best agreement with the experimental results of Gay et al. Q is obtained 
assuming a somewhat deeper groove (120 nm) as in the experiment (100 nm). This value is within the uncertainty of 
the actual milled depth using focused ion beam (FIB) fabrication. We have also presented a simplified model in which 
the surface wave is excited by a line dipole parallel to and just above a flat surface without groove. This approach 
permits the extraction of an analytical expression for the z component of the electric field along the interface. The 
agreement between this model and the full simulation is very good, showing that the transient "near-zone" regime 
does not depend on the precise shape of the groove. Indeed the details of groove depth and width only influences 
the amplitude and the phase of the generated wave. The overall form results from a line dipole source with a broad 
q spectrum interacting with a surface that supports a mode. Moreover, we have studied the influence on the wave 
structure of the propagative and the evanescent contributions. The propagative waves contribute importantly in 
the first few micrometers from the source, but their amplitude decay is faster and their wavelength is longer than 
the evanescent contribution. The wavelength of the propagative contribution decreases with the distance down to 
the excitation wavelength, whereas the effective wavelength of the evanescent contribution increases up to the SPP 
effective wavelength. Finally, we have studied a simplified model of the diffraction process, in which the reflection 
coefficient is replaced by a pole located at the SPP wave vector in q space. The scattered field can be then expressed 
in closed form. This "minimal model" correctly reproduces the SPP excitation and the wavelength evolution with 
distance. Such a semi-analytical model may be useful for the design and optimization of more elaborate photonic 
circuits whose behavior in large part will be controlled by surface waves. 
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APPENDIX 

The purpose of this appendix is to show location of the dipole, just above or below the metal/free-space interface 
is independent of the results obtained in section IIIII We begin with Eq. 

E z (x) = -r^ dqq[R(q) - R(oo) + i?(oo)l singx (12) 

-'"/••.I Jo ' 

The integrand diverges when q — > +00. However, the integral is defined for x > 0. We can write write Eq. (|12p as: 

Pa 

where 



E z (x) ~ 2 fc2 ) / dq qR(q) sin qx + R(oo) j dqq sinqa;, f (13) 



i?(oo) = ^4 and R = R(q) ~ R(oo) (14) 

1 + e(ujj 



so that from Eq.[U 



1 + eM VcMfcg-^+eHv^o^F 
The first part of Eq . [TBI converges because: 

qR{q) > -, with 7 G C. 

The second part is equal to the derivative of a Dirac function. Hence: 

E z (x) — ^—7-0 / dqqR(q) sinqx + aS'(x) (15) 
27rfc J 

The right term of the sum is only a term located at x — 0. For numerical integration, it is more convenient to use 
this last expression. 

If the dipole is located just under the interface: 

r+oo 

Ez( x ) — 7TJ2 / dqqT(q) sinqx, 

with: 



_ 2^Mfc 2 -'z 2 



Here. 



Hence: 



Ve(w)fcg - + e(w) ^Jk\ - q 2 



Tfo)— — >T(oo) = - 2 



{r-\-oo 
/ dqqf{q) sinqx 
Jo 

+T(oo) J dqq sin qxj 
dqqT(q) sinqx + j35'(x) (16) 



2?rfcg 



with 



T (q) = R(q) 



Hence the field diffracted by two dipoles located just above or just under the vaccum/metal interface differs only in 
a; = 0. 
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